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Abstract 

We construct analytic solutions describing black holes and black branes in asymptotically Lifshitz 
spacetimes with arbitrary dynamical exponent z and for arbitrary number of dimensions. The 
model considered consists of Einstein gravity with negative cosmological constant, a scalar, and N 
U{1) gauge fields with dilatonic-like couplings. We study the phase diagrams and thermodynamic 
instabilities of the solution, and find qualitative differences between the cases with \ <z<2, z = 2 
and z > 2. 



1 Introduction and summary of results 



One of the directions in which the AdS/CFT correspondence ^ has been extended in recent years 
is towards the construction of gravity models conjectured to be dual condensed matter systems with 
anisotropic scaling [21 [3l S]. It is expected that this direction of research, if successful in finding 
a dictionary translating between the gravitational degrees of freedom and field theory operators, 
would shed light onto the non-perturbative dynamics of non-relativistic models with this kind of 
scaling. 

In this work we will be interested in the development of the gravitational dual description of 
models exhibiting anisotropic scale invariance of the type 



where z is called the dynamic exponent. For z = 1, the scaling is isotropic; it corresponds to 
relativistic invariance. For generic values of z, the system is said to have Lifshitz scaling. The dual 
boundary field theory is not relativistic, but still allows for particle production. For the special case 
of z = 2 there is an extension of the Lifshitz scaling symmetry group to the Schrodinger group, in 
which particle number and special conformal transformations are conserved. See e.g. [H O [5l El [TJ [8] 
for some more references on this topic. 

In [1] it was proposed that gravity duals of field theories with Lifshitz scaling should have metric 
solutions that asymptote the form 



which is the generalization of anti-de Sitter spacetime (z = 1) to non-trivial dynamic exponent. 
A metric that locally looks like (|1.2p . we call a Lifshitz metric; it is invariant under the Lifshitz 
scaling if we let the radial coordinate scale as r — )• A~^r. The real parameter i represents the radius 
of AdS when z = 1, and we will refer to it as the radius of Lifshitz spacetime. 

In a holographic fashion, the dual theory would be formulated in a d-dimensional hypersurface 
located at infinite radial distance, the boundary. A boundary system on Mj x S'^~^ or Mt x M'^ with 
finite temperature and chemical potential can be realized from charged black holes or black branes 
respectively. In this paper, we will study the properties of asymptotically Lifshitz black holes and 
black branes by their own, and will not pay much attention to the holographic connections to non- 
relativistic field theory. We expect to address this issue in the future, and will have this application 
in mind throughout the paper. 

There is a large number of papers available in the literature discussing black holes in asymptoti- 
cally Lifshitz spacetimes, mostly using numerical methods. Some examples of analytic solutions, not 
necessarily with an Einstein gravity action, exist for fixed value of z [01 [TUl [TTl [T^ [T^ [HI [TBI [TU l [T7] . 
Numerical methods can be also employed to study a continuous range of z for both black holes and 
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black branes, as done for example in [THI [T9l [201 [2T1 [22l [23l [24] . 
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Most of the papers cited in the previous paragraph are based on an action containing a massive 
gauge field. In [25], an Einstein-Maxwell-scalar action with U{1) gauge invariance was considered, 
thereby providing an alternative realization of a gravity model supporting Lifshitz geometries. Both 
formulations have their own virtues, but also their disadvantages. The formulation with massive 
gauge fields can be embedded in supergravity models and string compactifications [5l[26 t [27 t [28l [29| 
[HKSl], but as mentioned above, the study of black holes is mostly numerical or for special values 
of the dynamical exponent z. This makes it less suitable or practical for applications in holography. 
On the other hand, the background in [25], including a dilaton-like scalar, is under better analytic 
control, and analytic charged black hole solutions can be found easily, as we demonstrate in this 
paper. The disadvantage of this model is that the dilaton is not constant and diverges on the 
boundary, which might be problematic for holography. It also means that the boundary theory is 
not exactly Lifshitz, but obeys some generalized scaling behavioi0. Presumably, a proper embedding 
of this model into string theory will shed more light on this issue. 

In this paper we will give analytic solution for black holes and black branes in any number of 
spacetime dimensions d > 3. This solution will have z > 1 as a continuous parameter and, in a 
sense, can be seen as a Reissner-Nordstrom version of asymptotically Lifshitz black holes. The 
matter content of the system we study, which is described by the action given later in equation 
(jl.3p . consists of > 1 abelian U{1) fields and one real scalar. It is an extension of the model 
considered in |25j , and is similar to the models studied in [30^ [3T] . Some of the gauge fields (how 
many depends on the topology of the black hole studied, as we will see below) and the scalar field 
are needed to support the Lifshitz spacetime at the boundary, and the remanent matter fields will 
contribute to Reissner-Nordstrom-like terms in the metric, which will have a clear signature in the 
thermodynamics of the system as a chemical potential term (in the grand-canonical ensemble). 

Before summarizing our results in the rest of this section, we find convenient to remind the reader 
of previous results found in asymptotically AdS (z = 1) spaces, both for the sake of comparing the 
solution to known cases and as a basis for our analysis. 

In [32] it was shown that, for black holes with spherical topology in asymptotically AdS space- 
times, there is a phase transition at a given temperature from a description in terms of ther- 
mal AdS (lower temperature) to a black hole setup. This Hawking-Page transition is due to a 
competing effect between the scale set by the volume of the spacetime and the scale set by the 
temperature. In [33] it was generalized to arbitrary dimension, and an explanation in terms of a 
confinement /deconfinement transition via a dual holographic field theory was given. The existence 
of unstable small black holes and the corresponding Hawking-Page transition, for the specific value 
z = 2, was predicted in [24] in a setup derived from string theory. We will find that the transition 
is present in the range 1 < z < 2 in our model. 

A U{1) gauge field was included in the setup in [34|, and the properties of charged black holes 
described by an Einstein-Mawell action with negative cosmological constant were computed. These 
solutions describe an asymptotically AdS Reissner-Nordstrom black hole with a horizon topology 

^We thank Marika Taylor and Kostas Skenderis for a discussion on this issue. 
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S at fixed time. Tlie existence of cliarged solutions gives rise to a rich pliase structure botli in 
the canonical and grand-canonical ensembles, which we describe below. 

The scale determined by the volume of the S"^~^ plays a crucial role in the construction of the 
phase diagrams, as in the case of the Hawking-Page transition commented above. A black brane 
solution (with a horizon with fixed-time topolog y W^-^) can be obtained by considering an infinite 
volume limit, which we will denote r/ — )• oo, where 77 is a dimensionless parameter to be introduced 
later. In this case the phase structure becomes trivial, with thermodynamics dominated by black 
holes for all temperatures. 

The results in j34j can be generalized to consider the abelian group U{1)^ without a significant 
change in the phase structure. Of course, the theories with M < N gauge fields can be recovered by 
setting N — M charge densities, pi, to zero. In figure [1] we outline a diagram showing the relations 
between theories with different number of gauge fields and horizon topology. 
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Figure 1: Diagram showing the relations between theories with one gauge field less, with boundary 
topology Mi X S'^^^ (BH, after black hole) and Mt x M'^"^ (BB, after black brane). The square 
represents asymptotically AdS spaces. The diagram continues indefinitely to the right and to the 
left up to U{1)^ , consisting on the Schwarzschild-AdS solution. The framed quantities indicate how 
many charges are there in the black hole. 



1.1 Summary of results 



In this paper we study the thermodynamic properties of charged asymptotically Lifshitz black 
holes. The solution we present is analytic for any value of the dynamic exponent z >! and for any 
number of spacetime dimensions d > 2. The model we consider is given by 
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(1.3) 



with electric fields in the radial direction turned on. Notice that this is a diffeomorphism-invariant 
action and we look for a metric solution that asymptotically approaches (jl.2p that breaks this sym- 
metry explicitely. In order to accommodate a Lifshitz spacetime in Einstein gravity, the presence 
of extra matter fields is required. The case with A'^ = 1 was studied for the first time in [25] . 
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where a BB solution was found. Strikingly, the solution for the metric in this case does not present 
the usual Reissner-Nordstrom properties that one finds in asymptotically flat/AdS solutions of the 
Einstein+Maxwell action. In concrete, the blackening function 6(r) defined in (|2.7p (defining the 
position of the horizon by means of b{rh) = 0) has just one non-negative root, and extremal 
black holes with finite entropy cannot be constructed. 

The gauge field is completely determined by the other fields present in the theory. In principle 
one would expect to have one free parameter associated to the gauge field, given by the constant 
of motion associated to A, since it appears in the action only through its derivatives. However, 
the requirement of having an asymptotically Lifshitz manifold (i.e., for z ^ 1) forces a relationship 
between this constant of motion and the magnitude of the scalar field 0. In a way, the role of the 
gauge field is to provide the appropriate potential to support an (asymptotically) Lifshitz spacetime, 
and the charge associated to the gauge field translates in the asymptotic properties of the manifold, 
and not in the horizon structure of the black hole. 

In the next section we present a generalization of this construction with an arbitrary number 
of U{1) fields. Remarkably, in the black brane case all the Aiy>2 gauge fields contribute to give a 
form of b{r) resembling that of the Reissner-Nordstrom solution. Indeed, the constants of motion 
associated to these extra gauge fields remain undetermined and are interpreted as charge densities 
Pi. In general, there is more than one root of 6(r) and extremal solutions with a finite horizon size 
do exist. As in the asymptotically AdS case, we can relate theories with a different number of U{1) 
fields by consistently turning off charges. In figure [2] we present a sketch of these relations. 

[7(1)^-1 C/(l)^ ?7(I)^+i 



A^- 1 




N 





Figure 2: Diagram showing the relations between theories with one gauge field less and the AdS 
limit for black brane solutions of our model. Circles represent asymptotically Lifshitz spacetimes 
whereas squares correspond to asymptotically AdS ones. The diagram continues indifenitely to the 
right and to the left up to U{1)^. The framed quantities indicate how many charges contribute to 
give the Reissner-Nordstrom-like factor in 6(r). 

Consider now the case with N = 2 gauge fields. As we just commented, it is possible to find 
a black brane solution with a rich horizon structure. One would be interested to know whether 
it is possible to find the correspondent solution but with a black hole, i.e., considering spherical 
symmetry in the form of the metri(|§. As we will show explicitely in the next section, such construc- 

^We are interested in this case and not in the hyperspherical one because it will describe finite size effects in an 
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tion leads to an algebraic equation whose solution fixes completely the form of both gauge fields 
in terms of the amplitude of the scalar field. The b function has a form reminiscent of that of the 
Schwarzschild-AdS case. This situation is completely analogous to the one found in |25j for the BB 
solution. In fact, the second gauge field also diverges at the boundary. Therefore, the role of this 
gauge field is to support the "sphericity" of the solution, and not to modify the horizon structure 
of the black hole. 

For the spherically symmetric case with more than 2 U{1) fields, the extra ^Ar>3 will again 
contribute in 6 to a new term resembling the charge-term in the Reissner-Nordstrom solution, thus 
modifying the horizon structure. As we did above, we present in figure [3] the relations between 
theories with different number of gauge fields when we turn off charges or take the AdS limit. 
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Figure 3: Diagram showing the relations between theories with one gauge field less and the AdS 
limit for black hole solutions. Circles represent asymptotically Lifshitz spacetimes whereas squares 
correspond to asymptotically AdS ones. The diagram continues indifenitely to the right and to the 
left up to U{1)^ . The framed quantities correspond to the number of charges contributing to give 
the RN-like factor of h . 

We have seen that the system described by the action (jl.Sp has a rich web of limiting cases. The 
question of whether we can recover black brane solutions by taking the appropriate t] ^ oo limit 
in the black hole solutions for asymptotically Lifshitz spacetimes has a positive answer. It turns 
out that this limit effectively decreases the number of U{1) fields in the theory by one, since now 
the inclusion of a gauge field supporting the "sphericity" of the solution is not required anymore. 
In figure m we present a partial diagram of the relations existing between the black hole and black 
brane solutions. 

1.1.1 Phase diagrams 

The study of the thermodynamic properties of the Lifshitz black hole solution leads to the phase 
diagrams showed in figure[5l both for the canonical and grand-canonical ensembles, where we sketch 
the results derived in the rest of the paper. 

In this figure we observe that the phase transitions depend crucially on the value of the dynamical 
eventual dual field theory. 
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Figure 4: Relations between the different specific cases arising from action (ll.Sp . See previous 
captions for symbolism. Tfie web of relations for the cases showed is incomplete for the sake of 
clarity, but can be completed with the previous diagrams. 



exponent z. For 1 < z < 2 we find a situation completely analogous to the one studied in [M] (which 
corresponds to the z = 1 AdS case). In the grand-canonical ensemble there is a line of first order 
phase transitions (the blue line on the top-left diagram), where the thermodynamically preferred 
solution is given by thermal Lifshitz spacetime at low values of the temperature and the chemical 
potential and by the black hole solution in the rest of the parameter space. As in the AdS case, 
the T = line is dominated by Lifshitz spacetime below a critical value of the chemical potentical 
<I>c (see equation (j4.58p ). and by extremal black holes for larger values. These extremal black holes 
have a non-vanishing entropy, and therefore are not expected to correspond to the true ground 
state of the theory. In the canonical ensemble, for low values of the charge, there is a first order 
phase transition between small and large black holes (for low and large temperatures respectively) 
which ends at a critical point (blue line and red point in the top-right diagram), above which the 
transition between small and large black holes is smooth. However, the solutions described in this 
paper present an electric instability, given by the shadowed region (see [35] for the AdS case). This 
instability implies that the addition of an infinitesimal charge to the black hole would contribute 
to a reduction of the electric potential. At zero charge a Hawking-Page transition between thermal 
Lifshitz spacetime and a black hole occurs, whereas at T = the solution is dominated by extremal 
black holes. 

When z = 2 the phase diagram in the grand-canonical ensemble remains qualitatively the same 
as in the 1 < z < 2 case. However, in the canonical ensemble the critical value of the charge for 
which there are no phase transitions above it is precisely given by Q = (see the middle-right 
diagram in figure [5]), and therefore only the Hawking-Page transition between thermal Lifshitz 
spacetimes and black holes setups remains. Furthermore, the electric instability also dissappears 
from the phase diagram at this value of the dynamic index. 

Last, for values of the dynamical exponent z > 2 the phase diagrams simplify even more. In 
the grand-canonical ensemble the line of first order phase transitions disappears, and black holes 
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z = 2 




4 



BH 



Q 



Q 



Tcril ThP 



BH 



z > 2 



Figure 5: Sketch of the phase diagrams obtained here for different values of the dynamic exponent 
z. See text for details. 



dominate the thermodynamics in almost all the parameter space. The exceptions occur at T = 0, 
where the extremal black hole dominates everywhere but exactly at the point $ = <I>C5 where 
the description is given by Lifshitz spacetime. Obviously, when <1> = the solution is also given 
by Lifshitz spacetime, since it is the only description available to us. In the canonical ensemble 
black holes dominate the phase diagram everywhere. These situations are shown in the bottom 
diagrams of figureEl where the black dots signal the points where the description is given by Lifshitz 
spacetime. 
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1.2 Outline 



The rest of the paper is organized in the following way. We start in section [2] by presenting the 
simple case in which the action (jl.3p has two U{1) gauge fields and the solution is of the black 
brane type. Once this is done we proceed to consider the case of the black hole with N >2 gauge 
fields, and analyze the limiting cases commented on in the previous summary. 

In section [3] we calculate the relevant thermodynamic quantities associated to the black holes 
found in section [2j These will be later used in section H] to unravel the phase diagram of the system 
for generic number of dimensions and dynamic index. 

We conclude in section [5] with some last comments and remarks. 

2 Setup and solution 

Consider the action (II. 3p . consisting of the usual Einstein gravity with a cosmological constant 
A < 0, a scalar, and U{1) gauge fields Ai coupled to the scalar. The equations of motion 
following from it are 

2yV 1 1 ^ f 1 "\ 

^^--J^9^.v--^d^cl>d,<^--Y.^e^''^\^^^^^-^^^^^^ = 0' (2-4) 

I)M(e^"'^^r) = 0> (2-5) 
N -I 

□0_^_A,e^.</'i^2 = 0. (2.6) 

i=l 

We are interested in finding an asymptotically Lifshitz spacetime. Furthermore, we will assume 
no dependence on the spatial directions, which will have the topology of M*^"^ in the black brane case 
or S"^~^ in the black hole case, and we restrict to the static case in which all the fields' dependence 
is along the radial direction of the asymptotically Lifshitz spacetime. We now choose an ansatz for 
the metric based on a single function bk{r), 

d.^ = - h,[r)^dt^ + r'dnl,_, , (2.7) 

with dri| the metric of a unit-radius S'^"^ if = 1 or the metric of M.'^^^ if /c = 0. ^ is the 
radius of the spacetime. One can also consider the case k = —1, based on a hyperbolic metric for 
di7|^_]^ but later on in this paper we derive that /c = — 1 is not acceptable unless z = 1. The 
function should asymptote 1 at large values of the radius to recover Lifshitz spacetime. The 
spatial components of the gauge fields are chosen to vanish to preserve rotational symmetry, and 
we will work in the gauge in which the radial component is also vanishing, therefore we will be 
concerned only with 74j^j(r) turned or|f| . Finally, for the scalar we consider it to be a function of 

■^To avoid confusion in our notation, let us emphasize that Ai^t{r) denotes the temporal component of the i— th 
gauge field A^,^, and not its time-derivative. 
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the radial coordinate only. 



2.1 Black brane with two f/(l) fields 

We consider first the warmup exercise in which we set k = (constant-r shces with Rj x M^^^ 
topology) and = 2 in the ansatz given before 

= ^b^) - + ^dxLi . (2.8) 

The Maxwell equations and the combination of Einstein equations — are solved by 

g0 = ^rV^id'Diz-i) ^ (2.9) 
{FiU = p.r-('^-^)e-^'^ (2.10) 

with i = 1, 2 and pi integration constants. These are related to the constants of motion associated 
to the gauge fields Ai, which enter in the action via the radial derivative, and therefore there are 
two conserved quantities 



5A'^ t 167rGrf+i ' 

which correspond to charge densities, as can be seen by calculating the total charge 



(2.11) 



with Vd-i = ^^""^J d'^~^x a dimensionless volume factor. 

The Einstein equation in any of the spatial directions, with the expressions (|2.9p - (j2.10p plugged 
in, gives a first-order differential equation for bo with solution 



{d-l){d + z-l) 2{d-l)^^d-z-l + y/2{d-l){z-l)X, ' 



The integration constant m will be related to the mass, as we discuss in the next section. With 
this solution at hand the rest of the Einstein equations, and the equation of motion for the scalar, 
become the algebraic equation 



4AV2(d-l)(z-l) = J] p2^-2id-i)-x.^2id-i)iz-i)^-x, _ _ ^2id-l){z-l)] . 

i=l 

(2.14) 

The r.h.s. of ()2.14p can be equal to the l.h.s. if we choose for the first gauge field 



A. = ->/2^ = -4A.^'f='-='^ ^ (2.15) 
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With these values, the contribution coming from the second gauge field has to vanish, which is the 
case if 

Let us comment what just happened. In order to satisfy the equations of motion we had to 
fix not only the coupling constant Ai, but also the charge of the first gauge field in terms of the 
scalar field amplitude /i. As discussed in the introduction, this gauge field is needed to support the 
structure of (asymptotically) Lifshitz spacetime. On the other hand, the second gauge field has 
a free charge, which will contribute to the thermodynamic analysis as a single chemical species. 
It is also the term responsible of having a bo function resembling that of RN black holes. With 
expressions ()2.15p - (l2.16p . the solution reads 

= ;3-5^-f27'«*npdC, , (2.17) 



A'l^j = rV2(rf + z - l){z - I) ^I'J ^ ^ (2.19) 

A'2,t = P2/i"^r2-'^-% (2.20) 
= /xrV2('i-i)(-i) ^ (2.21) 

where we have used 

2£2 ' V ^ 

to get the right asymptotics at infinity. 

It is now straightforward to check that in the uncharged limit, p2 ^ 0, one recovers the result in 
[25], whereas in the AdS limit, z — t- 1, the ^i^j field vanishej^. Without it, this solution is nothing 
but the AdS-RN solution considered in [3l]. These limiting cases are the ones outlined in figure [2j 

The solution presents a singularity at the origin r = 0, where curvature invariants diverge (except 
when z — )• 1 and p2 = 0, where the spacetime is AdS). However, the existence of an event horizon 
at a position r/j > cloaks it. The parameter m, which is related to the mass of the black brane, 
has to be positive definite, otherwise there will be no horizon and the singularity becomes a naked 
one. Following an argument on [36], we will impose the null energy condition T^^^^''^^ > with 

= (^/^, ^J-g^\Q) a null vector. From our solution it follows that T^ui^^" oc i'^{Rl - R\) = 
{d — l)(z — l)6o) so the null energy condition translates into z > \. This range of values of the 
dynamic index also ensures that we deal with real fields. 



^Actually, it vanishes provided that e'^ — ^ < 1. However, in the AdS hmit Ai — >■ — oo and the first U{1) field 
decouples from the system. 
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2.2 Generic case 



A question that arises after the analysis performed in the previous section is whether we could 
have obtained a black hole solution, i.e., a solution with k = 1. In this case equation (12.14p would 
have changed to (j2.24p below (with N = 2 plugged in). A solution does exist by fixing the charge 
of the second gauge field in a similar way as done before. The resulting black hole is similar to a 
Schwarzschild black hole in asymptotically Lifshitz spacetimes. 

In this section we will consider the black hole case with N > 2 gauge fields, which contains the 
case we just referred to. The metric is given by the ansatz (12. 7p with k = 1. As in the previous 
case, the Maxwell equations and a combination of the Einstein equations have as a solution for 
the scalar and the gauge field the expressions (j2.9p and (j2.10p . The first difference appears in the 
equation of motion for bi, which now reads (we keep an explicit factor k for later convenience) 

_ k{d-2) f 2Af ^-{d+z-l) ^ p2^-X,^o ^2(l-d)- V2(d-1)(.-1)A. 

''~d + z-3r^ {d-l){d + z-l) ^ 2{d - 1) d- z - 1 + ^2id - l){z - l)X, ' 

(2.23) 

With this, the algebraic equation equivalent to (|2.14p is 

AT 

= ^ p2^-2(d-i)-AV2(rf-i)(.-i)e-A,0o ^(rf _ 1);^. _ y^2{d - l){z - 1)] ^^C^^D 

i=l 



+2^/2{d -l)iz- 1) {k{d - l){d - 2)r-2 - 2A) . (2.24) 

As before, the first gauge field can be used to cancel the term proportional to A by choosing (|2.15p . 
and the contribution from the next N — 2 gauge fields can be cancelled by choosing 



^i = \hjZJ^ j = 2,... ,iV-l . (2.25) 



Finally, the A^-th gauge field can be used to cancel the term proportional to k if one fixes 



d — 1 \ z — 1 d + z — 3 

Notice that, once again, the charge of the first gauge field is fixed to support the Lifshitz spacetime. 
Additionaly, the A-th gauge field's charge is also fixed, in this case to support the existence of the 
S'^~^ topology for /c = 1. The hyperbolic case, k = —1, leads to imaginary charge densities (unless 
z = 1). Therefore, here and below, we will only consider the cases k = or k = 1. 
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Summarizing, the generic solution we found to (jl.3p is 

5^72" ~ £2 



n2 j„2 „22 

= ^^-^ft^dt^ + ^'dl^L-i, (2-27) 



/"^ ] + + y ,-2(c^+^~2) (2.28) 



J=2 



d-l 



yl'it = r V2(d + z - l)(z - 1) r'^+^"^ , (2.29) 



= p,/i"v2i^r2-'^-% (j = 2,... ,7V-1) (2.30) 



1 ,A/2A;(d-l)((i-2)(z-l) , <"-'^) 4 
yd + z — 3 

/ir^/2{53T)(iri) ^ (2.32) 



where again A = -{d + z - l){d + z - 2)/2l'^. 

We see that this solution depends on the parameters pj^ /i and m, which will correspond to the 
charge densities, the amplitude of the scalar field and the energy of the black hole, respectively. 
As commented several times already, two of the gauge fields have their charges fixed to support 
a spherical black hole in Lifshitz spacetime. The metric presents a horizon which, in general, has 
a near-horizon geometry given by the direct product of a 2-dimensional Rindler spacetime (the 
coordinates being r and t) and the spacetime given by d^k4-i-, which are spectator coordinates in 
this approximation. For a certain value of m and p the black hole becomes extremal. In this case 
the near- horizon geometry is given by AdS2 x S'^^^ or AdS2 x W^~^ depending on whether /c = 1, 0. 

We can take the AdS limit z — t- 1. In this case we recover once again the results in |34j . The 
fixed-charge gauge fields Ai^t and A^^t are set to zero (provided ^ < 1) and their coupling to the 
scalar field in the action goes to Ai^^r — t- — c«, decoupling them from the rest of the matter fields. 
The scalar field becomes constant as well. Figure [3] represents this limit and the reduction of the 
number of charges in consideration. 

We are also interested in recovering the solution characterized by A; = (this is, with flat 
topology of the constant-r slices) as an explicit limit of ()2.27p - (j2.32p . Following [3l], we introduce 
the dimensionless parameter r] and we scale r ^ Tjr. Given the form of the metric (j2.27p we must 
impose also t — )• ri~^t. We will focus on the neighborhood of a point in the S"^~^, considering 
just a fiat metric around it £'^ddf^_'^ — )• •q^'^dO,'^ Taking the scalar field to be scale-invariant, 
the appropriate scaling is fixed to be ^ — t- jy-\/2('^-i)(^-i)^. Now, the scaling on p determines 
completely the scaling of the Ai^t and A^r.^ gauge field charges. For the remaining pj charges we 
use the fact that these parameters actually correspond to charge densities, and therefore, to have 
scale-invariant charge, they have to transform under the scaling as pj — )■ 'q'^~^pj. All in all, the A'^ 
field strengths scale as 

Fi^Fi , ^ , Fn ^ v'^Fn (2.33) 
so in the large r/ — >■ 00 limit we should not consider the gauge field An, consistently with setting 
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/c = in the general solution. Furthermore, in this limit the metric becomes (|2.17p with a blackening 
function 69 given by (j2.18p (there N = 2). So we consistently reduced the black hole case to the 
black brane one, and are able to construct the web of relations depicted in figure HI 

An interesting special class of solutions (for k = 1 and = 3 for simplicity) are those that 
satisfy the mass-charge relation 



2 2(d - 2)2^2(1+.) - , 

= (ci- l)(d + . -3)3 • (^•^^) 
In this case, the black function reduces to 

5, = 1 + f - _^dS£=.-(^--)) ' . (2.35) 

\d + z-3r ^2{d-l)id + z-3) J 

This function is always positive and therefore there is no horizon. At r = there is instead a naked 
singularity. For z = 1 this is precisely what happens for the supersymmetric RN solutions. For 
z ^ 1, one may expect that the solution with (|2.34|) and (|2.35|) can be embedded as BPS solutions 
in some gauged supergravity action. These would need to be extensions of the present action (jl.3p . 
since in gauged supergravity with scalar fields, the scalar potential is not constant. 

Lastly, notice that even when the case d = 2 is not included in our analysis, it turns out that 
the solution also works in this case. For this value of the dimension, the solutions with A; = and 
k = 1 coincide locally. This is so because in this case there is only one spatial direction, and the 
difference between the two solutions is a global matter: whether the direction is compact or not. 



3 First law of thermodynamics 

We will discuss now the thermodynamic properties associated to the solution presented in ()2.27p - 
(j2.32p . For the sake of clarity we will restrict to 3 U{1) gauge fields with a black hole with spherical 
topology, i.e., just one non-trivial charge p2 = p. We will mark the difference between the black 
hole and black brane analysis keeping explicit terms of k, though. 



Temperature and entropy Unfornutately, it is not possible to obtain a general, analytical 
expression for the position of the horizon (given by the larger positive root of bk{rh) = 0) as a 
function of the three parameters m, p and p. However, we can still proceed to a thermodynamic 
study. 

Let us first notice that even when the parameter m is related to the mass of the black hole, 
as we will show in (I3.43p . it is not a fundamental parameter of the theory. These parameters 
are given by //, the amplitude of the scalar, p, corresponding to the charge (potential) in the 
canonical (grand-canonical) ensemble and the temperature T. Therefore we find it convenient to 
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express m = m{fi, p,T). However, as it is not possible either to obtain a closed expression for the 
temperature, we will use instead m = 171(^1, p,rh), given by 



m = r 



d+z-l 



(3.36) 



Notice that m is non-negative. Using this expression, the temperature as a function of the position 
of the horizon radius rh, p and p reads 



{d + z-l) + k 



-2(d+z-2) 



d + z-3rf^ 2{d-l) 



The entropy is given, as usual, by the Bekenstein-Hawking formula 

c _ ^d-l d-l 

The temperature (13.370 vanishes when the horizon radius satisfies 



p^ = 2id - l)W^^r- ( (d + . - 1) + k^^^ ] r2f , 



2 e2 



(3.37) 



(3.38) 



(3.39) 



'd + z-3r^^t^ 

where we have denoted with r^xt the position of the horizon at extremality, defined by the conditions 
b{rext) = b'{rext) = 0. Using the relation ()3.36p we can determine the value of tn^xt 



rriext 



k 



2{d-2fe r\ 



d + z-2pV V^d-1^22 



ext 



+ 



ext 



{d + z-^Yd + z-l d + z-3 {d-l){d + z-l) 



(3.40) 



Mass To evaluate the mass we should study the renormalized one point functions, where a coun- 
terterm must be added to regularize the expressions at infinity, since divergences will appear. 
However, here we wil take an alternative approach. We decide to evaluate the Komar integral 

1 



(3,41) 



with Kt = dt and substract the result from the thermal case (this is, the case with m = p = 
0, but such that the euclideanized time is periodic). We must match the normalizations of the 
Killing vectors between the black hole and thermal cases at Too to ensure they have the same norm 
expression at infinity. This is done by considering 



with Kq the Killing vector in the thermal setup. The result is 



M = Mt - Mq 



Vd- 



—ml-^-'{d-l) 



+1 



(3.42) 



(3.43) 
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where Vd-i is the volume of the unit S'^ ^ sphere. 

Expression (j3.43p will be useful in the grand-canonical case, where the charge is free to vary 
but the potential is fixed. In the canonical ensemble, however, one must fix the charge, which is 
proportional to p, and therefore the correct comparison scheme is to substract the result of the 
extremal black hole with the appropriate value for r^xt- It is not difficult to show that in this case 

AM = Mr - Me^i = (m - me.t)r^''{d - 1) . (3.44) 

Charges and chemical potential As commented above, there are three conserved charges 

= l6^^ ' ^^-''^ 
but two of them are completely specified in terms of the metric and the scalar 

= ,a r V2(rf + ^-l)(^-l)/^"V^ , (3.46) 
167rGrf+i 

The potentials associated to these charges in the thermodynamic relations come from the form 
of the the gauge fields potentials as functions of the radial coordinate 



S d + z-Y \ 



^3, = k ^^^^^^ 3^]) ^^^ v-^'^-^K^-i) (r^+^-^ - , (3.51) 

where we have fixed the integration constants such that the gauge fields vanish at the horizon, 
impliying that their norm-squared is non-singular there. The Ai and A^ fields diverge at the 
boundary, however they will not be of importance for a thermodynamic analysis, as we will see 
shortly. 

First law of thermodynamics It is now straightforward to check that the first law of thermo- 
dynamics hold^ 

dM = rd5 + $dg , (3.52) 



^From now on wc will denote Q2 = Q, since this is the only charge of importance in the thermodynamic relations. 
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where we define 



* = -*v(oo) = f^.i-'. (3.53) 



Specifically we have 

'dM\ ^ fdM 



^=\^] ■ (3-54) 



dSjQ \dQ 

In case we compare to the extremal case and not the thermal one, we have to consider the 
thermodynamic relation 

d(AM) = TdS + ($ - ^ext) dQ . (3.55) 
We can calculate now the heat capacity at constant charge, finding the result 

(dS\ dM/dru 

from where we conclude that the specific heat at constant charge is always positive and regular for 
z > 2. For 1 < z < 2 there is an instability for some values of the black hole parameters in the 
spherically symmetric case. We will comment further on this in the next section. 

Notice also that, for the first law of thermodynamics to be satisfied, the gauge fields Ai and 
(and their associated charges) are not needed. The fact that these two gauge fields do not 
seem to affect the thermodynamics may be related to having their charges completely determined 
by the scalar parameter /x and their diverging at the boundary, which would affect drastically a 
potential holographic interpretation. Therefore, it seems natural to assume that these fields, which 
are needed just to support the structure of the asymptotically Lifshitz spacetime solution, do not 
have a thermodynamic interpretation. 



4 Phase structure 

4.1 Grand-canonical ensemble 

Let us define the free energy from the thermodynamic relation 

W = M -TS -^Q , (4.57) 

where we have not included the contribution coming from ^\Q\ + ^sQs, i.e., as these two charges 
are fixed to support the asymptotic topology of the black hole solution, they must correspond to 
an ensemble in which those terms do not contribute to the free energy. In other words, we keep the 
charges Qi and Q-^ fixed, since otherwise these gauge fields would spoil the asymptotic topology 
of our solution. This is equivalent to keep the value of fixed in the phase structure analysis 
performed in the following. 
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Defining a critical potential given by 



2_,^2(d-l)(d-2)2 2(1-.) -./^ 



^ = k 

" {d + z-3)^ 



z-l 
d-1 



we can write parametric equations for the temperature T and the free energy W 
W 



.d+z-l 



z n-l-z 

An 



-z + {2- z) 



2{d-\) 



2\ i)2z-2 



(4.58) 

(4.59) 
(4.60) 



The z = 2 case Let us start the analysis of the parametric equations (j4.59p and (|4.60p by studying 
the case with z = 2. The temperature will be a bijective function of the radius of the horizon, and 
therefore every temperature is described by only one black hole. Clearly, there is a minimum value 
for T given by the setup with a vanishing radius of the horizon 



T, 



{rn=0,z=2) 



{d-l) 
8tt 



/iV^-i^($2 



$2) 



(4.61) 



The previous expression is negative when <I>^ > <I>^, and a quick look at equation (I4.60p shows 
that the temperature diverges as — )• oo. Therefore, being the relation between T and a 
bijective one, all positive values of the temperature are supported by a black hole. That this is the 
thermodynamically preferred solution is clear from the expresion for the free energy as a function 
of T and $ 

Vd-i f 47r 



W 



(4.62) 



On the contrary, for potentials less than the critical potential (in absolute value) the low temper- 
ature description has to be given by thermal Lifshitz spacetime. The line along which this phase 
transition occurs can be given in analytic form 



a/9;/ yad-i , / , , X 1/2 

$ = ^{d-2Y-An{d-i)T^ ~ - t] 



(4.63) 



where T, 



{z=2) 



{d - 2f/A'K{d - 1)1. 



The 1 < z < 2 case In this case, the competition between the two terms in the free energy (j4.59p 
depends on the sign of — This case is reminiscent of the z = 1 case studied in |34j . 

For values of the potential larger than the critical one the free energy is clearly negative and 
the black hole setup is favored in all the temperatures where such a description is valid. Now, in 
view of (j4.60p . it is clear that there are two terms: one positive and proportional to and one 
negative and proportional to There will be a positive value of the horizon radius at which 

the temperature vanishes, and above this value the temperature will be a monotonically increasing 
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T 

1.2 r 

1.0- 
0.8 - 
0.6- 
0.4- 
0.2- 

0.0 0.5 1.0 1.5 2.0 " 

Figure 6: Temperature as a function of the horizon radius for $ = l.l'&c (blue), ^ = 1.5<I>c (red) 
and ^ = 2^c (yellow). There is a minimum radius at which T = 0. An evaluation of the free energy 
shows that it is negative everywhere along the lines, specifically at T = 0. This plot was made for 
d = 6, z = 1.6 and /i = 2 in units where i = 1. 




function of the radius (see figure [6]). These two results imply that, for values of the potential 
greater than the critical one, there exists a black hole description of the system which, in turn, 
is thermodynamically preferred. Indeed, when T = at a finite radius of the horizon and the 
potential is not at its critical value the free energy is given by 

y,_ir^-- d-i 

SnGd+i d + z-3'' ' ^ ' 

showing that in this case the system at zero temperature is described by an extremal black hole, 
since W < when r/j > 0. 

When the potential ^ is less than the critical value the expression for the temperature (j4.60p is 
always positive and diverges when r/i — )• and ^ oo. When the horizon radius for a given po- 
tential is given by = 2z(d-i)^+~z^i) ~ *^^)' temperature approaches a minimum 
value 

1 2/2 

«22 / ^2 



_ d + z-l 
27r(2 — z) 



z){d + z-2,f 



-n^' — e^ ($2 -$2) 



9 (A IVW^ ^ y^c ^ ; ■ (4-65) 

2z{d — l)[d + z 

For any temperature larger than this there are two possible horizon radii. When we study the 
free energy, it is direct to see that W{Tmin) > 0. Two branches depart from this point, one 
corresponding to values of the horizon approaching r/j = and one approaching rh = oo, i.e., small 
and large black holes respectively. The branch corresponding to the small black holes gives always 
a positive free energy, whereas the branch associated to the large black holes corresponds to a 
negative free energy for cill va,lues of the tempercitiire T ^ 

(see figure [7|). Therefore, there is a 
first order phase transition between a thermal Lifshitz spacetime and a black hole when < 
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Figure 7: Temperature as a function of the horizon radius for $ = 0.2<^c (blue), $ = 0.5$c (red) 
and $ = 0.9<5c (yehow). There is a minimum temperature the black holes can describe. The dashed 
line marks the position at which the evaluation of the free energy along the lines changes sign, and 
it finishes just above the blue line, on the graph for $ = (not plotted). For small black holes the 
free energy is positive whereas for the larger ones it is negative. The values used to produce this 
plot were d = 6, z = 1.6 and ^ = 2 in £ = 1 units. 

The z > 2 case For values of the dynamic index larger than 2 and potentials less than the critical 
one (in absolute value), the free energy is negative for any value of the horizon radius, as can be 
seen directly from (|4.59p . Analyzing the expression for the temperature we observe that there are 
two positive terms multiplied by a positive power of r/j. Therefore, the temperature is a bijective 
function of the horizon radius, with vanishing temperature for zero horizon radius. Altogether this 
means that for low potentials and z > 2 the black-hole solution dominates the phase diagram, even 
at zero temperature. 

On the other hand, if <I>^ > the temperature will vanish at a finite horizon radius, in a similar 
way to the one reported on figure[6l Equation (j4.64p is still valid, showing that at low temperatures 
the thermodynamically favored description is given by the solution with a black hole. Indeed, this 
result holds true for any value of the temperature. 

The = case In all the previous cases we studied the phase diagrams for values of the 
potential above and below <I>c- When the potential is tuned to precisely its critical value the free 
energy can be expressed as a function of the temperature only 



At T = the preferred phase is the one described by the Lifshitz spacetime with finite potential. At 
any other temperature the black-hole description is the favored one. This signals the point T = 0, 




d+z-l 



(4.66) 
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^2 = $2 a special one, since its description is always given by Lifshitz spacetime, independently 
of the value of the index z. 

4.2 Canonical ensemble 

We now proceed to study the case in which we keep the charge Q (x p fixed, corresponding to the 
canonical ensemble. In this case the free energy is defined as 



F = AM-TS 



Vd- 



«— 1 — 2 r 



me^t{d-l)-zr''+'-' + k 



d+z-l I r ('^ ~ 2)^(2 — z) 2 d+2_3 



{d + z-3) 



(4.67) 



2{d-l){d + z-3)' 

with ruext given in expression (I3.40[) . As in the grand-canonical case, we should investigate sepa- 
rately the l<z<2, z = 2 and z > 2 cases. 

The 1 < z < 2 case Analyzing the temperature from p.37p for fixed charge, we observe that it 
can present an inflexion point when plotted against r/j. This happens at a position r^rit when the 
charge has the specific value pcrit 



{2-z){d-2fe 2 ^ 2z((i-l)(rf + z-l)^V'^-ir^^ 2(d+2-2) 

z{d + z-2){d + z-l) ' {d + z-3){2d + z-4) 



For values of the charge p < pcrit there is a region of temperatures described by three different 
black holes, with different horizon radii. If the charge is larger than the critical value then the 
relation between temperature and black-hole radius is in one-to-one correspondence. This situation 
is analogous to the one encountered for z = 1 [Ml , and can be seen in figure [5J In that figure 
we present also the results for the free energy, which show that the branch with dT/drh < 0, for 
charges below the critical one, is unstable, and a first order phase transition is present for these 
values of the physical parameters. At p = pcrit the kink in the free energy disappears and for larger 
values of the charge there is no phase transition. 

Another way to study this case is to consider the heat capacity at constant charge derived in 
()3.56p . The heat capacity associated to the unstable branch turns out to be negative, whereas the 
heat capacity for the two other branches, even beyond the kink in the free energy, is positive. In 
fact, one could overheat or undercool the system, keeping it in a metastable phase. In this case, the 
heat capacity grows as we enter further into the metastable region, eventually diverging. Figure [9] 
shows this behaviour. 

The z = 2 case When z = 2 something peculiar happens. Looking at the expressions for the 
critical radius and charge (j4.68p . we see that it occurs at the origin of spacetime, and that the critical 
value corresponds to the uncharged black hole. Any charge will be greater than pcrit and we do not 
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Figure 8: (Left) Temperature vs. horizon radius for p = 0.2pcrit (blue), p = pcnt (red), and 
p = 5pcrit (yellow). The parameters used are d = 4, z = 1.6, p = 2 and £ = 1. The position 
of the critical radius is marked with the dashed black line. (Right) Same color code for the free 
energy as a function of the temperature. The branch with dT/dr^ < in the left hand side plot 
corresponds to the unstable phase given by the cusp in the free energy plot. When T — )■ the free 
energy approaches from below. 

expect a phase transition. This can be seen by noticing that the term taking into account finite size 
effects in ()4.67p has disappeared. Therefore, for this particular value of the dynamic exponent, the 
free energy expression is equal to the one corresponding to the black brane solutiorj§. This means 
that, effectively, there are not two scales to compare (namely the horizon size and radius of the 
S'^^^) and therefore no phase transitions are present (notice however that the expression for the 
temperature still feels the finite size effects, so this is not quite the same as the study of the planar 
black brane). The heat capacity is positive for every value of the charge, signal of the absence of 
thermodynamic instabilities in the system for this value of the dynamical exponent z. 

The z > 2 case For z > 2 the term in the free energy proportional to k{z — 2) has reversed 
sign. One can show that the temperature is a bijective function of the horizon radius, which has 
a minimal value at which the black hole becomes extremal (the situation is completely analogous 
to the one shown in figure [6]) . The free energy is negative everywhere except at T = where it 
vanishes, and there is no non-trivial phase structure. Equivalently, the heat capacity at fixed charge 
is always positive. 

5 Final comments 

We have determined that for asymptotically Lifshitz spacetimes with a charged black hole in the 
center of the geometry, there is a phase structure that depends crucially on the dynamical exponent 

^Except for the factor Vd-i, which in this black hole case is finite, whereas in the black brane setup is infinity. 
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Figure 9: Heat capacity for black holes with charge less than, equal to and greater than the critical 
charge. The coloring is the same as in figure El We have split the plot in two graphs for clarity. 



z. 



For values \ < z < 2 the situation is reminiscent of that of the asymptotically AdS case, 
which corresponds to setting z = \. In the grand canonical ensemble there is a line of first order 
phase transitions. For large values of the temperature or the potential the thermodynamics are 
determined by the black hole configuration, whereas for low values of both parameters the ensemble 
is dominated by Lifshitz spacetime. When the temperature vanishes, for values of the potential 
larger than the critical one given in (|4.58p . the description is given in terms of extremal black holes 
with a nonzero entropy density. 

In the canonical ensemble there is also a line of first order phase transitions between two black 
hole setups, ending at a critical point given by the coordinates 

-2 



crit 



2 

Pcrit 



TT z{2d + z - 4)e 
2z{d-l){d + z - 



2-z 



z{d + z -2){d + z ~ 1) 



(2 



z{d + z- 2){d + z- 1) 



d+z-2 



(5.69) 



(5.70) 



(d + z-3)(2(f + z- 4)^22 

For values of the charge greater than pcrit there is no phase transition. For p = a Hawking- 
Page transition exists, with the thermodynamic ensemble dominated by Lifshitz spacetime at low 
temperatures and by the black hole solution for large temperatures. The exact point at which this 
phase transition takes place can be determined by calculating the non-trivial radius of the horizon 
at which the uncharged solution has vanishing free energy, and then plugging in equation (|3.37p . 
obtaining 



Thp 



d-1 



27r(2 - z)£ 



{2-z){d-2f 



21^/2 



(5.71) 



z(d + z-3)2 

Notice that the position of the critical point in the p — T plane depends crucially on z. Indeed, 
for the special case z = 2 one cannot find signs of a phase transition of any order at a finite 
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Figure 10: Free energy for the zero charge case with parameters (i = 4, /x = 2, £ = 1 and z = 1.5 
(blue), z = 1.6 (red), z = 1.75 (yellow), and z = 2 (dashed green). The kink is situated at T = Tmin 
and the curve with large gradient touches the F = axis at T = Thp- For larger values of z the 
free energy is negative for all values of the temperature except T = 0, where it vanishes. There is 
a branch with zero free energy from T = up to a finite value of the temperature. Then there is 
a kink at F = which gives rise to a branch with positive values of the free energy, until Tmj„ is 
reached. The right-most kink (not showed) gets closer to Tf/p as z is increased, until it coincides 
when z = 2. 




value of the charge, and only the Hawking-Page transition remains at a temperature given by 
Thp,z=2 = (d — 2)^/47r((i — 1)£. At precisely this value of the dynamic index, the temperatures 
Tmin and Thp coincide. To understand why, it is useful to plot the free energy for zero charge and 
various values of z. We do this in figure [TOl There we see that the value of the temperature at which 
the free energy presents a kink and is positive is given by Tmin, whereas this non-trivial branch 
crosses the horizontal axis at Thp- When we increase the value of the dynamical exponent z, the 
kink gets closer to the T-axis, and at precisely z = 2 it sits on top of it. If one increases further 
the value of the dynamical exponent, the free energy is negative for any value of the temperature, 
and therefore the Hawking-Page transition disappears, the uncharged solution being described by 
a black hole except at T = 0. 

For z > 2 there is a dramatic change in the phase diagrams, the phase transition disappears in 
both the canonical and grand-canonical ensembles, and the black hole description is the dominant 
one for any value of $,r (or p,T) except at the points {^c,T = 0), as given by equation (j4.66p . 
and (p = 0,T = 0). 

We have shown that the instabilities associated to thermal fluctuations, given by black holes with 
a negative heat capacity as given by equation (13.56p . correspond to thermodynamically unfavored 
phases. However, one should worry about electric instabilities as well. The isothermal susceptibility 
is given by 




(5.72) 



and has to be non-negative for the configuration to be stable. For 1 < z < 2, as was the case in [35] . 
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there is a region in the {Q,T) plane where such an electric instability is present. The line Qins{T) 
splitting the phase space into stable and unstable regions is given by the condition x{Qins,T) = oo, 
leading to' 



Pi 



'^^^^ {d-2f z{d + z-l) (2i:Tl{2 



d + 3 2 - z \ d + z-l 

(5.73) 

This curve exists only for 1 < z < 2, and encloses a region (for a given temperature, the values 
of the charge lower than Qins{T)) in which the solution presented here is unstable under electric 
perturbations. The critical point at which the line of first order phase transitions finishes is unstable. 
However, the temperature at which the Hawking-Page transition occurs is in the stable region, 
approaching it as z — )■ 2. 



Further directions 

We end with some comments for further study. As mentioned in the introduction, the Lifshitz back- 
ground supported by a dilatonic scalar suffers from divergencies at the boundary, which complicates 
a proper holographic formulation. Other issues related to Lifshitz spacetimes were discussed for 
example in [36^ I37j . It is expected that some of these problems will be resolved by studying more 
general models, for instance gauged supergravities with non-trivial scalar potentials that arise from 
string compactifications. Finding an embedding of our model into string theory is therefore worth 
investigating. 

The phases described in this paper should have an interpretation in terms of the dual boundary 
field theory. It would be interesting to identify and analyze the properties and phases of this field 
theory. For this, one needs to find an order parameter, which acquires different expectation values 
in the different phases. In the relativistic case, i.e. for z = 1, such an order parameter is given by 
the Wilson loop. 

It is also interesting to perform a thorough investigation of the critical exponents at the phase 
transition and thermodynamic instabilities, as was done in the AdS case in [38]. The critical expo- 
nents may be sensitive to the value of the dynamical exponent z, as might the properties of the 
universality classes. Another line of research goes along the work performed in [39], where transport 
coefficients were studied in asymptotically Lifshitz spacetimes. 

Finally, another extension of our model is to add bulk fermions in the probe approximation. 
This allows us to study properties of condensed matter systems with fermions that obey Lifshitz 
scaling. A first step in this direction is under investigation |40j . 
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